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Abstract. Any method for estimating the ensemble average of arbitrary operator 
(observables or not, including the density matrix) relates the quantity of interest to a 
complete set of observables, i.e. a quorum. This corresponds to an expansion on an 
irreducible set of operators in the Liouville space. We give two general characterizations 
of these sets. All the known unbiased reconstruction techniques, i.e. "quantum 
tomographies" , can be described in this framework. New operatorial resolutions are 
■ given that can be used to implement novel reconstruction schemes. 

\o ' 
o 
o 
o 



03 



X 



1. Introduction 

In order to characterize a quantum system, one can measure an observable, or a set of 
observables, on repeated identical preparations of the system. A question immediately 



arises: is this set a quorum, i.e. is it sufficient to give a complete quantum information 
on the system? In other words, is it possible to estimate the expectation values of any 



system operator? This issue is of a great practical interest for fundamental experiments 
in quantum measurement theory, as well as for a potential application as a quantum 
standard in the new field of quantum information 

Different estimation techniques have been proposed tailored to different systems, 
such as the radiation field @, trapped ions and molecular vibrational states H, spin 
systems [|J, and a unified approach is desirable. For set of operators that exhibit a group 
symmetry, a general theory, the so called "group tomography", has been established 
H . However, not all estimation techniques can be described within a group-theoretical 
scheme, and the purpose of this paper is to give a more general framework. As we 
will see, all the known quantum estimation techniques can be embodied in the present 
approach. In addition, the formalism here presented is of help for the derivation of new 
operatorial resolutions that result into novel estimation techniques. 

The paper is structured as follows. In section we establish the general conditions 
for a quorum of observables. Sections |3], £|, and || are devoted to examples of the 
estimation technique, for the harmonic oscillator, general spin systems, and the free 
particle respectively Section ^| closes the paper with a summary. We added [Appendix 
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A to reformulate the theory in the familiar Dirac formalism, and [Appendix B| to give a 
constructive algorithm to derive tomographic basis for finite-dimensional Hilbert space. 



2. Quantum estimation 

The indirect (tomographic) reconstruction f| of an operator A is possible when there 
exists a resolution of the form 



A 



/ dx Tr 


AB\x) 


Jx 





C(x) 



(1) 



where x is a (possibly multidimensional) parameter living on a (discrete or continuous) 
manifold X . The only hypothesis in ([]]) is the existence of the trace. The operators C(x) 
are functions of the quorum of observables measured for the reconstruction, whereas the 
operators B(x) form the dual basis of the set C(x). The term 



£[A](x) = Tr AB\x) C(x 



(2) 



represents the quantum estimator for the operator A. Indeed, the expectation value of 
A, namely the quantity of interest, is given by the ensemble average 



(A) = Tr 



Ag 



dx Tr 



x 



A B\z 



Tr 



C{x)g 



dx(S[A](x)) 



(3) 



x 



where g is the density matrix of the quantum system under investigation. The averaged 



estimator in Eq. (0) is the product of two terms: The quantity Tr C(x)g 



depends only 



on the quantum state, and it is related to the probability distribution of the measurement 



outcomes, whereas the term Tr 



AB\x) 



depends only on the quantity to be measured. 
In particular, the tomography of the quantum state of a system corresponds to writing 
Eq. (|1]) for the operators A = \k)(n\, {\ri}} being a given Hilbert space basis. For a 
given system, the existence of a set of operators C(x), together with its dual basis B(x) 
allows universal quantum estimation, i. e. the reconstruction of any operator. 

We now give two characterizations of the sets B(x) and C{x) that are necessary 
and sufficient conditions for writing Eq. (|1|). 

Condition 1: bi-orthogonality 

Let us consider a complete orthonormal set \n) (n = 0, 1, • • •). Formula (0) is 
equivalent to the bi-orthogonality condition 



x 



dx{q\B\x)\p) (m\C(x)\l) = 6 mp 6 lq 



(4) 



where 5ij is the Kronecker delta. Eq. 
to a continuous basis. 



can be straightforwardly generalized 



| Among the existing estimation techniques for quantum systems, the so-called quantum homodyne 
tomography of a single radiation mode has received much attention in the literature |s) . In this paper, 
the term "tomography" is collectively used to denote any kind of state-reconstruction technique. 
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Condition 2: completeness 

If the set of operators C(x) is irreducible, namely if any operator can be written 
as a linear combination of the C(x) as 



A 



dx a(x) C(x) 



x 



then Eq. (pi) is also equivalent to the trace condition 



Tr 



B\x)C{y) = S(x,y) 



(5) 



(6) 



where S(x,y) is a reproducing kernel for the set B(x), namely it is a function 
or a tempered distribution which satisfies 



dx B(x) S(x,y) = B(y) 



x 



One can easily see that an analogous identity holds for the set of C(x) 



dx C(x) S(x,y) = C(y) . 



(7) 



(8) 



A' 



The proofs are straightforward. The irreducibility condition on the operators C(x) is 
essential for the equivalence of (|l|) and @. A simple counterexample is provided by the 
set of projectors P(x) = \x)(x\ over the eigenstates of a selfadjoint operator X. In fact, 
Eq. (P) is satisfied by P(x). However, since they do not form an irreducible set, it is 
not possible to express a generic operator as 6 ^ f x dx (x\0\x) \x)(x\. 
If either the set B(x) or the set C(x) satisfy the additional trace condition 

B\y) B(x) 

~C\y) C{x) 



Tr 



Tr 



S(x,y) 
S(x,y) 



(9) 
(10) 



then we have C(x) = B(x) (notice that neither B(x) nor C(x) need to be unitary). In 
this case, Eq. ([!]) may be rewritten as 



A 



dx Tr 



x 



AC\x) C(x) 



In abstract terms a certain number of observables Q x constitute a quorum when 
there are functions f x {Qx) = C(x) such that C(x) form an irreducible set Notice that 
if a set of observables Q x constitutes a quorum, than the set of projectors \q) xx (q\ over 
their eigenvectors provides a quorum too, with the measure dx in Eq. (JJ) including 
the measure dq. Of course it is of interest to connect a quorum of observables to a 
resolution of the form (|]), since only in this case there can be a feasible reconstruction 
scheme. If a resolution formula is written in terms of a set of selfadjoint operators, the 
set itself constitutes the desired quorum. However, in general a quorum of observables 
is functionally connected to the corresponding resolution formula. If the operators 
C(x) are unitary, then they can always be considered as exponential of a set of 
selfadjoint operators, say Q x . The quantity Tr C(x)g is thus connected with the 
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moment generating function of the set Q x , and hence to the probability density p(q; x) 
of the measurement outcomes, which play the role of the Radon transform the quantum 
tomography of the harmonic oscillator M. Here, the basic resolution formula involves 
the set of displacement operators D(a) = exp[aa* — aa], which may be viewed as 
exponential of the field-quadrature operators x<$, = \(ae~ 1 ^ + a)e 1 ^). In general, the 
operators C(x) can be any function (neither self-adjoint nor unitary) of observables 
and, even more generally, they may be connected to POVMs rather than observables. 
In Sections |37| and ^]4| we will see examples of this situation. 



3. Quantum estimation for harmonic system 

The harmonic oscillator (HO) model provides a detailed description of several systems 
of interest in quantum mechanics, as the vibrational states of molecules, the motion 
of an ion in a Paul trap, and a single mode radiation field. Different proposals have 
been suggested in order to reconstruct the quantum state of a harmonic system. As we 
will see in the following they can be summarized using the framework of the previous 
section, which is also useful for devising novel estimation techniques. 



3.1. Quantum homodyne tomography 

Perhaps the most famous quantum estimation technique is given by the so-called 
quantum homodyne tomography @. Homodyne tomography applies to a single-mode 
radiation field as well as to the vibrational state of a molecule or a trapped ion, 
and consists of a set of repeated measurements of the quadrature operator = 
^(ae~ l( ^ '+ a) e 1 ^) at different values of the reference phase 0|| For homodyne tomography, 
the relevant operatorial resolution is provided by the set of (irreducible) displacement 
operators D(a) = exp (aa) — aa) ,a G C. In fact, for the displacements Eqs. (||) and 
(fLOl) hold since 

Ti[D(a)b^((3)} = 7rc^(a - (3) , (12) 

and Eq. ( |TTD reduces to the Glauber formula 

d 2 a 



A = / Tr 

7T 



AD* (a) D{a) . (13) 



Changing to polar variables a = (i/2)ke 1 ^, Eq. ( [IB] ) becomes || 

A = T— Tr(i e ik «*) e"^ , (14) 



71 4 
" J— oo ^ 

which shows explicitly the dependence on the quorum q^. After taking the ensemble 
average of both members, evaluating this trace over the set of eigenvectors of q^, one 
obtains 

A = r # r°dqp{q;<t>)R[A]{q;<l>) , (15) 

Jo * J-oc 

§ For the vibrational tomography the quadrature operator is a time-evolved position or momentum. 
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where p{q;4>) = 4>(q\Q\q) $ is the probability distribution of quadratures outcomes. The 
tomographic kernel for the operator A is given by 

R[A}{q-<P) = Ti[AK{q-q^)}, (16) 

where the integral kernel K(z) is 



K(z) 



lim -Re- — 

e^o+ 2 (z + is) 2 



;i7) 



P denoting the Cauchy principal value. 

Using our condition one can see that the Glauber formula can be generalized to 

d 2 a 



A 



Tr 



C 



AF 1 D(a)F 2 F^D^^F^ 1 



where F% and F 2 are two generic invertible operators. By choosing F± = F 2 = S(Q with 



S(() = exp 



i (C V 2 - <V) 



(eC 



the squeezing operator, we arrive to a different tomographic resolution 



(A)= - 
Jo * 



+00 



dq P{(q; <p) Tr A K[q 



(19) 



(20) 



in terms of the probability distribution of the generalized squeezed quadrature operators 



= i [(fie* + z/e-^y + (juT** + ue**)a] , 



(21) 
(22) 



where /i = cosh |£| and z/ = sinh |£| exp(2i arg[£]). Such an estimation technique has 
been investigated in detail in Ref. K 



3.2. Phase-space estimation techniques 

A different estimation technique based on the generalized Glauber formula ( 18|) may be 
obtained by putting F\ = I, the identity operator, and F 2 = (— ) flta , the parity operator. 
In this case one gets 

d 2 a 



A 



Tr 



7T 



AD\a){- 



l D(a) 



(23) 



and therefore changing variable to a = 2(3 and using the relation (-) a a D(2f3) 

d 2 [3 



(A) 



Tr 



c n 



A4D\(3){-) ata D{/3) Tr D((3) g b\(3) {-) aU . (24) 



Eq. ( p4j ) says that it is possible to estimate a HO operator A by repeated measurement 
of the parity operator on displaced versions of the state under investigation. An 
approximated implementation of this technique for a single mode radiation field has 
been suggested 0, [12|] through the measurement of the photon number probability on 
signals displaced by means of a beam splitter. A similar schemes has been used for 
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the experimental determination of the motional quantum state of a trapped atom [13[ . 



The advantage of ( g4j ) compared to the approximated methods is in the possibility of 
directly obtaining the kernel if [A] (a) for any operator A for which the trace exists. For 
instance, the reconstruction of the density matrix in the Fock representation may be 
obtained by averaging the kernel 

K[\n)(n + al\\}{a) = (n + d\4 D j (a){-f a D{a)\n) (25) 
= 4 (-)<+* exp{-2M 2 } J^yy (2a)^(4|a| 2 ) , (26) 



without the need of artificial cut-off in the Fock space [13 . 



3.3. Nonlinear phase tomography 

Let us now consider the set of selfadjoint operators given by 

0) = V^(V) M0) Hl>) ^ , E [0, 2vr) , (27) 

where //(</>) is the canonical London' POVM describing the ideal measurement of the 
HO phase and V(ijj) describes a nonlinear phase-shift, namely 

M0) =|e*)<e*| 

= exp [i (a f a)V] , (28) 

where |e i?i ) = Xm e I 71 ) * s t ne Susskind-Glogower vector. For a single-mode radiation 
field the action of V(ip) corresponds to a nonlinear Kerr interaction, whereas for a 



trapped ion it could be obtained by laser excitation of vibronic levels [14" . 

One may argue that the measurement of the operators B(ip,4>) provides a complete 
characterization of the state under investigation. This is indeed the case with B(ip,<j)) 
as a self-dual basis, as it can be easily proved using the bi-orthogonality condition (f|) 

" 27r f 27T dS dib 

I 77-77- (m\B(ifj, <f>)\l) (q\B(ip,<j))\p) = 5 p 2_ q 2 !m 2_ l 2 5 p _ q ^ m = 5 mp 5i q . (29) 
Jo ^ ^ 7r 

The kernel for the operator A is obtained from equation (|3|) 

K[iK0,V>) = Tr [i%>,V)] • (30) 
In particular, the kernel K[P n d]((j), ip) for the matrix elements is given by 

K[P nd )((f),ip) = (n + d\B(<f),ip)\n} = exp [#(d 2 + 2nd) + i(j>d] . (31) 
Notice that, for diagonal matrix elements, Eq. (31) needs a regularization procedure. 



In fact, the kernel for the projector P nn is given by 

K £ [P nn ] (0, i>) = exp [i2ipne + i<pe] , (32) 

with the limit e — > that should be taken after the average over the probability density 
p(4>, ip) = Tr[g B((p, tp)}. A similar procedure should be employed for the reconstruction 
of any operator which is a function of the number operator a^a only. 
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The nonlinear "phase tomography" presented in this section represents a novel 
resolution formula for harmonic oscillator operators. The corresponding reconstruction 
technique is based on an ideal phase measurement of the Kerr-displaced state. Hence, 
in this scheme the quorum is a POVM. 



3-4- A nonunitary resolution formula 



In this section we present a resolution formula in term of nonunitary operators and show 
how its implementation would correspond to a generalized measurement i.e. a POVM. 
The operators R n (4>) with neZ and G [0, 2n) are defined as follows 



Rn{<t>) 



e+ e ^ 



n > 
n < 



4(0 



e -ia>a^ e -n n<Q 



(33) 



where e_ = |n + l)(n\, e + = J2 n \ n )( n + 1| are the so-called lowering and raising 
operators for the harmonic oscillator. The operators R n {4>) are not unitary, however 
Rn(4>) plays the role of dual basis for R n {4>). This can be easily seen using condition 
(|), for instance we have 



Tr 



n,k>0 



g n e ia>a(f 



Tr \e~ ia1a ^ e k _ 



(34) 



and similarly for the other cases. As a consequence a generic operator A may be written 

as 



(35) 



For estimating HO operators one needs a recipe to obtain the expectation value of R n (4>). 
Using a resolution of identity in terms of phase vectors |e^) 



I 



2tt 



one evaluates the traces as 



w d± p -iq4> l(k _ 



Tr 



gR g (ip) 



fir d4 



j: n ^e^ + ^ (4>\g\4> + ^) q<0 



(36) 



(37) 



Eq. |37| implies that the knowledge of (R n (4>)) is equivalent to that of the density matrix 
in the phase representation. 



4. Quantum estimation in spin systems 

The recently born spin tomography ^ |J allows to reconstruct the quantum state of a 
spin s system. It employs measurements of the spin component in different directions, 
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i. e. it uses as quorum the set of operators S-n, where S is the spin operator and n a unit 
vector. Various different quorums may be constructed by exploiting different directions. 

The easiest choice is to consider all possible directions. The procedure to derive 
the tomographic formulas for this quorum is analogous to the one employed in Sect. 
3TTJ for homodyne tomography. The reconstruction formula for spin tomography for 
the estimation of an arbitrary operator A, using the measurement outcomes m of the 
component of the spin in all directions n = (cosy? sin sin (p sin cos $), is 

( A ) = E / % PK*0 ^[i](m,n) , (38) 

m=-s ^ n 

where p(m, n) is the probability of obtaining the eigenvalue m when measuring S ■ n, 
R[A](m, n) is the tomographic kernel for the operator A, and Q is the unit sphere. 
In this case the operators C of Eq. ([[]) are given by the set of projectors over the 
eigenstates \m,n) of the operators S ■ n for all directions n. Notice that it is a set of 
irreducible operators in the system Hilbert space 7i. In order to find the dual basis 
B, one must consider the unitary operators obtained by exponentiating the quorum, 
i.e. D(i[),n) = exp(iipS ■ ft), which satisfy the bi-orthogonality condition (|4]). In fact, 
D(ip,n) constitutes a unitary irreducible representation of the group SU(2), and the bi- 
orthogonality condition is just the orthogonality relations between the matrix elements 
of the group representation |J, i.e. 

J dg b jr {g)D\ k {g) = ^5 jk 5 tr , (39) 

where D is a unitary irreducible representation of dimension d, dg is the group Haar 
invariant measure, and V = j R dg. For SU(2), with the 2s + 1 dimension unitary 
irreducible representation D(ip, n), Haar's invariant measure is sin 2 ^ sin $d$d(pdi/j, and 



^ = j^zp Thus, the bi-orthogonality condition is 

J dnj^ # sin 2 t ( 3 \ e ^\r)(t\e-^\k) = 6 jk 5 tr , (40) 
Using Eq. (flCf) and condition 1, it is immediate to write the spin tomography identity 



A = — dn dib sin — Tr 

4vr 2 In Vn ^ 2 



2s + 1 f ^ f 2n , ; _ : _ 2 tfj 



AD*(i>,n) D(i>,n) . (41) 



Notice the strict analogy between Eq. fl4T|) and Glauber's formula (|T3|), analogy which 
derives from the group symmetry that underlies both homodyne tomography (Weyl- 
Heisenberg group) and spin tomography (SU(2) group). In fact, both these tomographies 
may be derived in the domain of Group Tomography ||. In order to obtain the 
reconstruction formula (^), one only has to take the expectation value of both members 
of Eq. (^T|) and to evaluate the expectation value trace on the eigenstates \m, n) of S ■ n. 
Thus, the explicit form of the tomographic kernel is obtained as 

R[A](m,n) = [ d^j sin 2 £ Tr \A e ~^")] . (42) 

7T Jq 2 L J 
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As already anticipated, there are other possible quorums for spin tomography. For 
example, for spin s = | systems, a self-dual basis for the operator space is given by 
the Pauli matrices and the identity {^o~ x , ■^d'y, ^^z, ^1}- Hence, we would expect 
to find a tomographic identity where the above operators constitute a quorum. In fact, 
from the property a a -a a > = I6 aa > (a, a' = x, y, z), it is immediate to see that both the bi- 
orthogonality relation ([|) and the trace condition hold. The following reconstruction 
formula for an arbitrary 2x2 matrix A derives 

i 

2 

(A) = Yl p( m '^>) Tr 

m=-\ ot=x,y,z 

In the case of generic s spin system, Weigert has shown J| that it is possible, by 
choosing (2s + l) 2 arbitrary directions for n, to obtain (in almost all cases) a quorum of 
projectors \s, nf) (s, nf\ (J = 1, • • • , (2s + l) 2 ), where \s, ftj) is the eigenstate pertaining 
to the maximum eigenvalue s of S ■ fij. The search for the co-basis is done numerically, 
but since we have a discrete quorum, it is possible to use the procedure derived from 
Gram-Schmidt orthogonalization method, which is presented in |Appendix B . 



m H — Tr 
2 



A 



(43) 



5. Quantum state of a free particle 

Can we infer the state of a moving packet from position measurement in time? The 
answer is positive as can be rigorously proved using our condition (|J). The probability 
density of the position of a free particle at the time r is obtained from the selfadjoint 
operator 

R(x,t) = e~ ip2T \x)(x\ e ip " T , (44) 

where \x) are eigenstates of the position operator and p the momentum operator. We 
suppose for simplicity a particle with unit mass and use normalized unit h/2 = 1, such 
that the free Hamiltonian is given by Hp = p 2 ■ The dual basis is constituted by the set 
of operator R(x, r) themselves, as follows from Eq. and from the following relations 
( Ij)? 3 — a i b, c, d denote position eigenstates) 





dx dr (a\R(x,r)\b) (c\R(x,r)\d) 

dxdr e -ir(b 2 -a*+<?-d?) ( a \x)(x\b) (c\x)(x\d) 

dx dT e - iT ( b2 - a2 + c2 - d2 ) e ix(a-b+c-d) 

= 5{a - c) 5{b - d) . (45) 
Therefore, a generic free particle state can be written as 

Q= / dx dr p(x,r) R(x,t) (46) 
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p(x, t) = Tr[g R(x, r)] being the probability density of the particle to be at position x 
at time r. Eq. ( fj6|) can be generalized to particle moving in arbitrary potential, for 
details see Ref. fl5| . 



6. Summary 

A unified approach to quantum state estimation has been presented. Some general 
conditions that guarantee that a set of observables is sufficient for the estimation were 
given. A constructive algorithm to derive dual basis, and therefore quantum estimators, 
is suggested for the finite dimensional case. Our framework allows to describe all known 
quantum tomographies (i.e. unbiased state estimation procedures) for single mode 
radiation field, trapped ions, vibrational states of molecules and spin systems. Moreover, 
using our characterizations, some new tomographic resolutions were derived. These may 
be used to implement novel estimation techniques. 



Appendix A. Liouville space formulation 

In this section we reformulate the general scheme of quantum estimation, given in Sect. 
0, by making use of the properties of the linear space of operators acting on the system 
Hilbert space 7i, i.e. the Liouville space C(7i). This may results in a more convenient 
visualization of of our conditions 1 and 2, since customary Dirac notation can be used. 
The Liouville space is built starting from the set of operators of Hilbert-Schmidt type, 
which is itself a Hilbert space, with the scalar product 

(A\B) = Tr \A ] B] . (A.l) 

Using Dirac formalism on the operator Hilbert space, we find that kets are the operators 
of £(H), i.e. O = \0), while bras are obtained with the substitution Tr[ 1 O^] = (0\ 1 . 
The space of operators may be extended by considering non normalizable vectors in 
C(Ti). Using Dirac notation, we can rewrite all the formulas introduced so far. The 
generalized tomography formula (p]) corresponds to the expansion of the vector \A) on 
the non-orthogonal basis \C(x)), that will be referred to as "spanning set", and on its 
dual (B(x)\ (x G X), i.e. 

\A) = [ dx\C(x))(B(x)\A) . (A.2) 

In this framework, conditions 1 and 2 represent respectively the identity resolution for 
the spanning set and the definition of the dual vectors {B(x) \ of the set \C(x)). In fact, 
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we can readily rewrite Eqs. (PS) and 



dx \C(x))(B(x)\ 



x 



(A.3) 

(B(x)\C(y)) = 6(x,y). (A.4) 

The linear structure of the Liouville space provides a necessary and sufficient 
condition to verify that a set of operators is a quorum. A set of vectors \C(x)), such 
that (B(x)\C(y)) = S(x,y) Vx, is a spanning set (with (B(x)\ as dual set) iff the only 
operator \0) G C{Ti) that is orthogonal to all \C(x))s is the null operator, i.e. iff the 
two equivalent conditions 



(6\C(x)) = Tr 6 ] C{ 







[B(x)\0) = Tr B\x)6 







(A.5) 



(for any x G X) imply that O = 0. 

In addition to the reconstruction of operators acting on the system Hilbert space 
H,, one can extend the formalism also to the reconstruction of super-operators acting 
on the system operator space. A typical example is the Liouvillian super-operator that 
evolves the system density operator into a density operator. In fact, by introducing two 
resolutions of the identity in the operator space, one can express any super-operator L 
in terms of its "matrix elements" (B(x)\L\C(y)) on a basis for £(H), i.e. 



L 



A' 



dx / dy\C(x)){B(x)\L\C(y)){B(y)\. 



(A.6) 



x 



By taking a basis {\n)} of 7i, Eq. ( |A.6|) rewrites as 



L k ± = (k\L \m){l\ \p) 



mp 



dx / dyTi C(x)\p)(k\ Tr B(y)\m)(l\ 
x Jx 



Tr 



B(x)L[C(y)} 



(A.7) 



Appendix B. Orthogonalization procedure 

Here we give an algorithmic procedure, usable in the case of finite quorums, to 
construct the set of dual operators (B n \ of the quorum \C n ). Using the Gram-Schmidt 
orthogonalization method, one obtains a basis \yk) from a complete set of vectors |C/-) 
(assume for simplicity that all |C%) are non-zero and that in {|Cfc)} there are no couples 
of proportional vectors): 

\yo) = i |G>> 

|| Eq. (A. 2) is readily obtained from Eq. (|]) by introducing the identity super-operator, defined as 
1[A] = A for any A £ £(H). In fact, by using the basis {|n)} for H, we see that 

(fc|i \m)(l\ \p) = S m kSi p . 



(B.1) 
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where N = || \C ) || and N k =\\ \C k ) - J2j=o \yj){Vj\ C k) II- Now, by writing the identity- 
resolution 

1 = £ \y k ){y k \ = l -^(y \ + £ i ( |C fc > - £ l%X%l^) ) <V*I, (B.2) 

fc=0 k=l k \ j=0 J 

and using repeatedly Eq. ( p.l|) (expressing \yf) of Eq. (p.2j ) in terms of the \C n )s and 
reorganizing the terms) we can find the dual set (B n \ as 



(B.3) 

Unfortunately no such a general procedure appears to exist for the case of a continuous 
spanning set. Many cases,though, satisfy the conditions @ and (0), and hence we can 

write (B(x)\ = (\C(x)) X 



Acknowledgments 

This work has been partially supported by INFM through project PAIS-1999-TWIN. 



References 

[1] H. K. Lo et al, Introduction to Quantum Computation and Information, World Scientific (1998). 

[2] G. M. D'Ariano, P. Kumar, and M. F. Sacchi, Phys. Rev. A 61, 013806 (2000) 

[3] T. Opatrny and D. -G. Welsch, Prog. Opt. vol. XXXIX, 63 (1999). 

[4] S. Weigert, Phys. Rev. Lett. 84, 802 (2000). 

[5] G. M. D'Ariano, Phys. Lett. A 268 151 (2000). 

[6] K. Vogel, H. Risken, Phys. Rev. A 40, 2847 (1989). 

[7] G. M. D'Ariano, C. Macchiavello, M. G. A. Paris, Phys. Rev. A 50, 4298 (1994); G. M. D'Ariano, 
U. Leonhardt, M. Paul, Phys. Rev. A 52 R1801 (1995). 

[8] G. M. D'Ariano, "Measuring quantum states", in Quantum Optics and the Spectroscopy of Solids, 
ed. by T. Hakioglu and A.S. Shumovsky, Kluwer Academic Publishers (1997), p. 175. 

[9] F. D. Murnaghan, The theory of group representation, pg. 216, Johns Hopkins Press (1938). 
[10] G. M. D'Ariano, S. Mancini, V. I. Manko, and P. Tombesi, Q. Opt. 8, 1017 (1996) 
[11] K. Banaszek and K. Wodkievicz, Phys. Rev. Lett. 76, 4344 (1996). 
[12] T. Opatrny and D. -G. Welsch, Phys. Rev. A, 55, 1462 (1997). 
[13] D. Leibfried et al, Phys. Rev. Lett. 77, 4281 (1996). 
[14] S. Wallentowitz, W. Vogel, P. L. Knight, Phys. Rev. A 59, 531 (1999). 
[15] Th. Richter, A. Wiinsche, Acta Phys. Slov. 46, 487 (1996). 



